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We investigate the self-organization of a binary mixture of similar sized rods and dipolar soft spheres by means of Monte-Carlo 
simulations. We model the interparticle interactions by employing anisotropic Gay-Berne, dipolar and soft-sphere interactions. In 
the limit of vanishing magnetic moments we obtain a variety of fully miscible liquid crystalline phases including nematic, smectic 
and lamellar phases. For the magnetic mixture, we find that the liquid crystalline matrix supports the formation of orientationally 
ordered ferromagnetic chains. Depending on the relative size of the species the chains align parallel or perpendicular to the 
director of the rods forming uniaxial or biaxial nematic, smectic and lamellar phases. As an exemplary external perturbation we 
apply a homogeneous magnetic field causing uniaxial or biaxial ordering to an otherwise isotropic state. 


1 Introduction 

In the last decades significant progress has been made in 
understanding the phase behavior of mixtures of rods and 
spheres. The great interest in these systems was coined by the 
pioneering work of Asakura and Oosawa in 1954^ in which 
they predicted that an effective attractive force arises between 
two large particles in a suspension of considerable smaller de¬ 
pletion agents. As a consequence, even the simplest mixtures 
composed of purely repulsive rods and spheres exhibit a rich 
phase behavioi^2. In the subsequent period, phase and reen¬ 
trant transitions as well as multiphase equilibria between gas, 
liquid, liquid crystalline (LC) and crystal phases have been ob¬ 
served W investigations based on experiments®, theoretical 
studies®® and computer simulationsEMHI. 

Correspondingly, even more complex behavior is expected 
for suspended particles with internal degrees of freedom such 
as magnetic particles. As predicted in a famous work by 
Brochard and de Gennes the coupling between the ferrocol- 
loids and the liquid crystal molecules i ncreas es the sensitiv¬ 
ity to external magnetic field significantl>M^ ^, a ph enomenon 
with many pot ential applications for magneto or electro- 
optical devicesIEini. So far, suspensions of ferrocolloids in a 
nematic LC matrix composed of much smaller LC molecules 
(also referred to as low molecular mass LCs) have been stud¬ 
ied experimentally and theoretically (see, e.g.. Refs. rfTM 
and references therein). 

At the present time, increasing attention is being paid to 
the study of mixtures of spherical magnetic and rodlike col- 
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loids where the spheres are of similar size as compared to the 
width of the rods. Recent experimental examples include mag¬ 
netic nanoparticle suspensions®, attempts® to produce sus¬ 
pensions of magnetic nanoparticles in dichroic pigment parti¬ 
cles® and percolating carbon nanotubes®. 

In the present study, we consider such a colloidal rod-sphere 
mixture with additional permanent dipole moments for the 
spheres. Dipoles exhibit a long-ranged dipole-dipole interac¬ 
tion with the head-to-tail alignment being the most energet¬ 
ically favorable spatial configuration. In fact, dipolar parti¬ 
cles alone show even in the absence of an external (electric 
or) magnetic field complex micro structure formation i nclud- 
ing chain-like®, network-like®, and ring-like structures®®. 
Prime examples of self-assembled dipolar structures fueled by 
external fields include chain formation in constant fields®, 
layer formation in rotating fields®, and the structure forma¬ 
tion in triaxial fields® with a wide rang e of applications in 
magnetorheology® and biomedicineEES^ 

The aim of the present study is three-fold: firstly, to set-up a 
tractable model for rod-(nondipolar) sphere mixtures and ex¬ 
plore the corresponding LC order, secondly, to provide infor¬ 
mation to which extend both types of ordering, i.e. ferromag¬ 
netic and LC ordering, can be combined in mixtures involving 
dipolar interactions, and thirdly, to illuminate how far the col¬ 
lective behavior can be manipulated by coupling an external 
magnetic field to the dipolar spheres. 

To our knowledge, no theoretical investigations have been 
carried out on the particle level for such a magnetic hybrid 
system and thus, a microscopic understanding, is so far miss¬ 
ing. As a first step, we have recently presented a Monte-Carlo 
(MC) simulation study® targeting two specific size ratios in 
the absence of any perturbation. In the present paper, we aim 
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to explore in more detail the behavior of the underlying ref¬ 
erence system (involving non-magnetic spheres), as well as 
other relevant size ratios. Moreover, we investigate the impact 
of the simplest external perturbation, that is, an homogeneous 
external field. 

The remainder of the paper is organized as follows. In Sec¬ 
tion 2 we formulate the theoretical model and give a brief out¬ 
line of the simulation details. Section 3 presents results for the 
rod-sphere mixture. We begin with the non-magnetic system 
in Section 3.1 followed by a discussion of the dipolar system 
in Section 3.2. The infiuence of an external magnetic field is 
presented in Section 3.3. We conclude in Section 4. 


For the strength anisotropy parameter in Eq. ^ we useEl 

e(u,-,Uj,f;j) = eo [ei(u,-,u;)]'' [e 2 (u,',u;,r,;)]^ (4) 

with Sq, ei(ui,Uj) = [l — X^(u,'^ being the strength 
parameters from the original formulation of the overlap model 
by Berne and Pechukas^S (where ji, v are adjustable expo¬ 
nents). The latter parameter in Eq. Q, e 2 (u;,Uy,f/y), has been 
introduced later by Gay and Bemel^to adjust the well depth 
ratio for the side-by-side (e^) to end-to-end (e^) configura¬ 
tion by introducing the parameter x' = + 

^y^), that is. 


2 Model 

Our model fiuid consists of a mixture of dipolar soft spheres 
with embedded, permanent dipole moments and uniaxial rods. 
The total interaction energy for the system may be decom¬ 
posed into contributions stemming from each species and a 
contribution accounting for interactions between rods (r) and 
spheres (s), that is, 

iVs A^s A/r A/r iVs A/r 

i=lj^i «=l7=l 

where is the number of spheres and Nj^ is the number of 
rods. 

For the interaction between rods we use a single-site poten¬ 
tial model suggested by Gay and BemeEl which considers an 
orientation-dependent range parameter between pairs of parti¬ 
cles. 

The latter reduces in the case of two identical particles of 
length I and width Oq tcP 


e2{ui,uj,rij) =1 


^ r (r,7-u,- + r,7-uy)^ 
2 [ I+Z'uru; 

1 -X'urUj 


( 5 ) 


For the pair potential between two dipolar soft spheres 
(DSS) of diameter cr^ we use 


where u^s is the (truncated and shifted, cf.l^ soft-sphere (SS) 
repulsion given by 

^ss ) = 4^0 I ‘^112 
\^ij\ 

and t/dd accounts for the dipole-dipole potential stemming 
from the permanent dipole moments m/, that is. 


Udd{mi,mj,rij) = -3-( 8 ) 
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Finally, for the interaction between pairs of rods and spheres 
we consider a Gay-Berne potential of the same form of Eq. 

To this end, we have to determine the range parameter depend¬ 
ing only on the orientation of rod 7 , uj, and the normalized 
(2) vector r^/ connecting to the centre-of-mass of sphere /, that 
is El] 


where X = ~ is the director along 

the principal axis of particle i and rij is the connecting vector 
between the centre of masses of particles i and j. With the 
range parameter a defined in Eq. 0 we can now introduce 
a modified Gay-Beme (GB) potential of Lennard-Jones form, 
that is, El 


UUui,Uj,rij) = 4e{ui,Uj,rij) 


<7o 


|r,7| -a(u;,Uj,f,7) + o'o 
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(To 
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( 3 ) 


(Trs(u 7 ,f,y) = G^q^ [l -/'(f.-y • Uj)^] ^ (9) 

where Gq = ^(ob + (Ts) and x!' — “ (Tq )/(/^ -f G^). Sim¬ 

ilarly, the strength anisotropy for pairs of rods and spheres 
becomes 


ers(u;,r,7) = £0 [l-x"\rij-Ujf]^ (10) 

where the rod-sphere well-depth anisotropy is given by x'^' = 
1 — (e^/e^)^/^. For our calculations we choose an established 
set of parameters for which a large number of simulation stud¬ 
ies have been carried out, that is, I/gq = 3, £s/^e = 5, /i =2 
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and V = 1 . In the present study the diameter of the spheres 
remains as an adjustable parameter with a* = g^/go. For the 
numerical computations of the potentials we use cutoff radii 
of = 4(7 o for the rod-rod interactions, = 2.5Gs for the 
short-ranged soft-sphere repulsions, and rf = 2{lGs) for the 
rod-sphere interactions. The long-ranged dipolar interactions 
are treated W the Ewald summation method with conducting 
boundaries!^. 

We have examined systems for a variety of particle com¬ 
positions (va = Nsi/N where a = s,r for spheres and rods, 
respectively), volume fractions ( 0 a = where t/a is a 

particles’ volume), and total number densities p* = Ng^/V. 
We have performed Monte-Carlo (MC) simulationsES in the 
canonical ensemble to examine the morphologies of these sys¬ 
tems. Equilibration requires of the order of 2x10^ cycles and 
a further 5x10^-1x10^ cycles is being used for the calculation 
of ensemble averages of quantities of interest. A MC cycle 
consists of N trial attempts (moves, orientations and moves- 
orientations) for a randomly chosen particle. 

The orientational ordering of the phases has been quantified 
with the aid of order parameters. More specifically, we mea¬ 
sure the orientational ordering via the nematic order parameter 
5*^^^ obtained through diagonalizing the ordering tensoi® 

Qlc = I [3 (^f )b (KHc - ^c] /2N, (11) 

1=1 

where b,c = v,y,z (cartesian components) and Kf is the carte¬ 
sian component of the symmetry axis of rods or dipolar 
spheres (along the direction of mp. For each species, the 
eigenvector associated with the largest eigenvalue 5+ of the 
ordering tensor is considered with respect to the director lis or 
fir. The other two eigenvalues, and S-, fulfill the inequality 
> So > S-. The biaxiality of the phase can be calculated 
via the order parametei® 

5= ^^(nr-ns)^-I) (12) 


which yields B = 1 (parallel directors) and 5 = — 0.5 (perpen¬ 
dicular directors) for an uniaxial and a biaxial phase, respec¬ 
tively. The polarity of the magnetic phase is measured via the 
first rank polar order parameter, that is. 



An analysis of translational order ing has been performed 
by means of correlation functionsS^ffil We have calculated 
(i) the longitudinal correlation function 



Lw^(n|-ko'-"a|) \ 

AV2P {N, - 1) j 


(14) 


of the projection of the intermolecular vector parallel 
to the macroscopic principal director ha where Ay 2 = 
71 — {vij • lia)^^ Ar|| is the volume of a cylindrical shell with 

thickness Ary =0.05, (ii) the dipole-dipole correlation func¬ 
tion perpendicular to the director fia, that is, 

(a) / Lw g ('-X - V rfj - (rp- - fla)") cos Gjj \ 

\ L'^7 s{r±-yjrfj - {nj ■ fia)^) / 

(15) 

where cos Otj = % • zj and z/ is the pricipal axis (given by the 
dipole vector). Finally, to analyze the structure with respect to 
the dipole moment of a DSS particle we have calculated the 
two-dimensional correlation function 


{q,r±) = 

I Zi^j 5 (r± - d (r,, - |r,j • m,|) 

\ AVp{N,-l) 

where AV = n + Ar^)^ — Arp and m, is the dipolar 
unit vector of particle i. 

3 Results and Discussion 

3.1 Binary mixtures of rods and soft spheres (GB-SS) 

In this section we study the self-organization of a binary mix¬ 
ture of Gay-Beme rods and repulsive soft spheres (GB-SS 
mixture) for different sphere diameters ranging from cr* = 1.0 
to ( 7 * = 2.0. The GB-SS mixtures are considered as reference 
systems since the particles’ self-organization is essential for 
the examination of more complex binary mixtures of rods and 
dipolar soft spheres (GB-DSS). In particular, the configura¬ 
tions of the GB-SS system are considered as initial configu¬ 
rations for the GB-DSS system. We have focused on systems 
with Xr = 0.8 and = 0.9 for which fully miscible phases 
are obtained. This is not a trivial finding since the amounts 
of spheres that can be supported by LC phases depends cru¬ 
cially on the chemical affinity of the species® Hence, binary 
mixtures of (GB) rods and spherical particles which attract 
each other, such as Lennard-Jones spheres, exhibit demixing 
transitions® or microphase separation occurs in hard-core in¬ 
teracting rod-sphere mixtures!^. Depending on the size of 
the spheres these values correspond to volume fraction ra¬ 
tios (rods to spheres) 2.4 < 0r/0s < 12 which is of the order 
of magnitude that can be achieved in real colloidal suspen- 
sions“ 
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Fig. 1 (a) Tentative state diagram of a GB-SS mixture with (a) 

CTs* = 1 and Vr = 0.8, (b) = 1.5 and Vr = 0.8, (c) G^ = 1.7 and 

Vr = 0.8 and (d) G* = 2 and = 0.9. The points on the diagram 
indicate the pairs (r*,p*) for which the actual simulations were 
performed. The solid and dashed lines indicate state transformations 
of the GB-DSS (magnetic) system. Abbreviations: Fully miscible 
isotropic (I), uniaxial nematic (Nu), uniaxial smectic (SmB) and 
lamellar (L) phase. 


3.1.1 GB-SS mixtures with a* = 1.0. Initially, we exam¬ 
ine GB-SS mixtures with cr* = 1.0. We have performed simu¬ 
lations for systems consisting of = 720 particles for Xr = 0.8. 
Larger systems of = 2000 particles have also been studied 
(for selected state points) to account for finite size effects. Ei¬ 
ther cooling series starting from a high temperature isotropic 
liquid or a nematic phase (prepared by melting a crystal state) 
are performed, at specific total number density p*. In some 
cases, heating series from ordered phases were performed as 
well to check for the reversibility of the phase behaviour. 

A calculated temperature-density (r*,p*) diagram is pre¬ 
sented in Fig.[^ where we use different symbols for the three 
distinct phases. Specifically, fully miscible phases of different 
order are exhibited: (i) an isotropic (I), (ii) an uniaxial nematic 
(Nu) and (iii) a smectic (SmB) phase. 


In the I phase both species lack long range orientational or¬ 
der. A characteristic snapshot of the isotropic phase is shown 
in Fig. ii- The orientational disorder of the rod species 
has been confirmed by calculating the order parameter 
A representative order parameter vs temperature diagram for 
p* = 0.4 is given in Fig. [^. By performing a cooling se¬ 
ries starting from an I-state the system first transforms into an 
uniaxial Nu state where the rods possess a long range orien¬ 
tational order as can be seen from the increase of the order 
parameter (see also the snapshots in Fig. |^). In the Nu state 



Fig. 2 Representative snapshots of a GB-SS mixture with G* = 1 
and Vr = 0.8 in various states: (a) isotropic state at 
[(r*,p*) = (1.8,0.40)], (b) uniaxial nematic (Nu) state at 
[(r*,p*) = (1.2,0.40)] and (c) smectic B (SmB) at 
[(r*,p*) = (0.8,0.40)]. The direction of fir is also shown. For 
clarity, the rod species have been removed from the simulation box 
in the bottom row. 


both rods and spheres are homogeneously distributed and the 
phase has uniaxial symmetry with respect to the nematic di¬ 
rector. The pair correlation function (ry) confirms that 
both species are distributed homogeneously along fir. 






Fig. 3 Order parameter as a function of temperature of a GB-SS 
mixture with (a) CTs* = 1.0 at p* = 0.40 and (b) G* = 2.0 at 
p* =0.338. Pair correlation function for the rod (r) and 
spherical (s) species with (c) G* = 1.0 at p* = 0.40 and (d) 

CJ* = 2.0 at p* = 0.338 for various T*. 


By further decreasing the temperature, a state with transla¬ 
tional order is obtained in which the rods are organized in peri- 
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odic layers (see Fig[^). This behavior is also confirmed by the 


:). The rods possess 


strong modulations in (r||) (see Fig ' 
local hexagonal order within the layers similar to monodis- 
persed (GB) rod systems®. Interestingly, in this SmB phase 
there is a tendency of the spherical particles to organize into 
linear arrangements along the director of the phase as it can be 
seen from inspecting the snapshots shown in Fig|^. This is 
also reflected by the behavior of the function (ry) which 


displays modulations within the SmB phase with a periodic¬ 
ity of approximately one Oq (see Fig. |^). This implies that 
the spheres are incorporated within the layers of rods rather 
than positioning between successive layers of rods through 
microphase separation. We discuss this issue in more detail 
in the next subsection. 


3.1.2 GB-SS mixtures with 1.5 < cr* < 2.0. We now turn 
to study GB-SS mixtures with larger diameter of a* > 1.5. 
Specifically, we have performed simulations for systems of 
N = 720 particles with cr* = 1.5 and cr* = 1.7 for = 0.8 
and of = 1251 particles with a* = 2.0 for = 0.9. Tem¬ 
perature vs density (r*,p*) diagrams are shown in Fig. [^-d. 
A common feature with the previous system (a* = 1.0) is that 
fully miscible I and Nu phases are obtained. Both I and Nu 
lack long range positional order. Typical snapshots of these 
phases are shown in Fig |^-b. An example that monitors the 
order parameter as a function of the temperature within the Nu 
state is given in Fig. for a system with cr* =2. 



Fig. 4 Representative snapshots of a GB-SS mixture with a* = 2, 
for Xr = 0.9 in various states: (a) isotropic state at 
[(r*,p*) = (2.4,0.338)], (b) uniaxial nematic (Nu) state at 
[(r*,p*) = (1.8,0.338)] and (c) Lamellar (L) state at 
[(r*,p*) = (0.8,0.338)]. For further details see Fig. 2. 


Remarkably, the layered phase that is obtained for lower 
temperatures is not a common smectic phase. The trans¬ 
lational ordered phase consists of alternating rod-rich and 


sphere-rich layers indicating spontaneous microphase sepa¬ 
ration (see Fig. [^). Particularly noteworthy is that both 
(^11) (^11) are modulated in the smec¬ 

tic phase revealing translational ordering along the director 
fir. It should be noted that there is also considerable inter- 
digitation especially for the rodlike particles. This type of 
smectic phase is termed lamellar (L) and has already been ob¬ 
served both experimentallyl^ and theoreticallyEI^II] in partic¬ 
ular, i mpl ementations of MC simulations® and Onsager type 
theory^ for purely repulsive (hard-core) rod-sphere mixtures 
have shown that the lamellar state is thermodynamically stable 
for cr* = 1.0. Interestingly, we do not observe lamellar orga¬ 
nization for the GB-SS mixture with cr* = 1.0. The reason for 
this is because, in our model, the side-by-side rod-sphere con¬ 
figuration is energetically favourable. Therefore, three spheres 
can fit into one layer spacing of rods without destroying the or¬ 
der of the SmB state. It should be noted that the side-by-side 
configuration is not favoured for rod-sphere mixtures interact¬ 
ing through hard pair potentials (such as hard spherocylinders 
and hard spheres®). In this case, in order to reduce the ex¬ 
cluded volume of spherocyfinder-sphere pairs the species mi¬ 
crophase separate and form a lamellar phase®. For systems 
with ( 7 * > 1.5 a side-by-side rod-sphere configuration is also 
the most energetically favourable one; in this case, though, 
this configuration destroys the smectic order of the rods since 
the spheres do not fit into one layer spacing of the rods. In 
order to preserve a layered structure at lower temperatures the 
system microphase separates into a lamellar phase with a high 
degree of interdigitation for the rod species. 


3.2 Binary mixtures of rods and dipolar soft spheres 
(GB-DSS) 

In this section we examine the role of the dipolar interactions 
on the structure of binary GB-DSS mixtures. The configura¬ 
tions of the non-magnetic reference systems (see section 3.1) 
are used as initial configurations for the GB-DSS systems. We 
have considered spheres with a relatively large central perma¬ 
nent point dipole moment of st rength m* = m/ SqCJ^ = 3. 
Inspired by experimental values®^^® the coupling parame¬ 
ter A = m^/ksTG^ takes values greater than 3.5. 

3.2.1 GB-DSS mixture with a* = 1.0. Initially we 
present a brief summary of the behavior of a GB-DSS mixture 
where the diameter of the spheres is equal to the width of the 
rods. A detailed analysis is presented elsewhere®. The topol¬ 
ogy of the (r* — p*) diagram is similar to the respective GB- 
SS system that is shown in Fig[^; nevertheless, the morphol¬ 
ogy of the DSS particles within the phases is completely dif¬ 
ferent. The solid and dashed lines in Fig[T^ indicate I-uniaxial 
nematic Nu and Nu-uniaxial smectic SmB transformations, re¬ 
spectively. An important finding is that the I-state is destabi- 
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lized in favor of the LC phases in the absence of any external 
stimuli. In the I-state, the DSS particles self-assemble into 
isotropic networks of wormlike chains. Remarkably, in the 
Nu state the ferromagnetic chains spontaneously align along 
the LC director fir. Overall, the ferromagnetic chains are ran¬ 
domly arranged “up” or “down” and the system does not ex¬ 
hibit a net magnetization. Finally, within the smectic state the 
rod particles form a SmB phase with well defined layers; the 
ferromagnetic chains penetrate and pass through these layers. 
This result is in accordance with the observed configurations 
of the corresponding GB-SS mixtures in which the SS parti¬ 
cles fit into the smectic layers (see Sec. 3.1.1). 

3.2.2 GB-DSS mixture with a* = 1.5 From here on we 
examine the role of the rod-sphere size ratio on the self¬ 
organization of GB-DSS mixtures. We address basic ques¬ 
tions regarding (i) the morphology of the DSS particles inside 
the LC matrix and (ii) how the DSS ordering influences the 
LC matrix and vice versa. As we have seen in Sec. 3.1, the 
translational order changes considerably by tuning the rod- 
sphere ratio already for the case of the (non-magnetic) GB- 
SS mixture. First, we increase the size of the dipolar spheres 
to (7* = 1.5. A temperature-density (T* — p*) diagram has 
been calculated for binary mixtures of A = 720 particles and 
Vr = 0.8 (see Fig[^). This corresponds to volume fraction ra¬ 
tios (rods to spheres) 0r/0s — 3.54 which is of the same order 
of magnitude that can be achieved in real colloidal suspen- 
sions^S]. Isotropic (I), uniaxial nematic Nu and highly inter- 
digitated smectic (Sm) phases are obtained. 
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Fig. 5 (a) State diagram of a GB-DSS mixture with a* = 1.5 
(xr = 0.8), involving isotropic (I), uniaxial nematic (Nu) and 
uniaxial interdigitated smectic (Sm). The symbols indicate state 
points where the actual simulations were performed, (b) Global 
order parameters as functions of T* at p* =0.387. 


In the I state both species do not possess long range orien¬ 
tational and positional order. The DSS particles self-assemble 
into ferromagnetic chains. A representative simulation snap¬ 
shot is shown in Fig [^. A simila r stru cture is also found in 
monodispersed dilute DSS systems^SEll gy reducing the tem¬ 
perature starting from the isotropic phase the system under¬ 
goes an I-Nu phase transformation. A significant increment 


of the order parameter 5*^^^ of both species is clearly demon¬ 
strated in a (5*^^^ — r*) diagram for p* = 0.387 that is shown 
in Fig 1^. Remarkably, the ferromagnetic chains are spon¬ 
taneously unwrapped within the nematic phase forming fer¬ 
romagnetic chains that are on average parallel to the director 
hr. Interestingly, there is a notable delay concerning the mag¬ 
nitude of the order parameter in comparison to 5*^^^ that 
jumps to 0.4 at the I-Nu transformation. This signifies a lower 
alignment of the ferromagnetic chains by the LC matrix in 
comparison to the GB-DSS mixture with cr* = 1.0 (cf. Fig. lb 
in Ref. [34]). The corresponding value of B 1.0 means that 
the directors of the species are, on average, parallel to each 
other rendering the phase uniaxial. A characteristic snapshot 
of the Nu phase is shown in Fig[^. A signature of the arrange¬ 
ment of the DSS particles are the periodic arcs that appear in 
the anisotropic correlation function (parallel 

to the dipole in comparison to the perpendicular direction) (see 
Fig|^). 



Fig. 6 Representative simulation snapshots of a GB-DSS mixture 
with a* = 1.5 for Xr = 0.8 in various states: (a) isotropic state at 
[(r*,p*) = (2.8,0.387)] (b) uniaxial nematic (Nu) state at 
[(r*,p*) = (1.6,0.387)], and (c) highly interdigitated smectic (Sm) 
state at [(r*,p*) = (0.8,0.387)]. For clarity, only DSS particles are 
shown in the bottom row. The directions hr and h^ are indicated by 
thick lines. 


Interestingly, the nematic phase is significantly enhanced in 
favor of the isotropic phase in comparison to the system with¬ 
out dipolar interactions (see Fig.[^ and Fig. [^). Therefore, 
there is an interplay in which the LC rods induce orientational 
order to ferromagnetic chains and vice versa. Even though the 
chains are polar, the phase does not show spontaneous magne¬ 
tization since the polar chains are arranged into an antiparallel 
manner. This has been confirmed by the polar order parameter 
Pi that takes a value near zero. The presence of polar domains 
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beyond a single chain is also excluded, since the {rj_) 
function takes large positive values near the origin and van¬ 
ishes already at small distances (see Fig|7j)). It is imp ortant 
to note that monodispersed systems of dipolar spheres®!!!, 
in the absence of external fields, do not exhibit spontaneous 
orientational order for the densities considered here. The ho¬ 
mogeneous positional distribution of rods along the director 
is confirmed by the structureless function (ry) (see Fig 

^). The light modulations of (ry) in the Nu state with 
periodicity of approximately one molecular diameter shown 
in Fig.[^ occurs due to the correlations of spheres that belong 
to the same chain; furthermore, the absence of a peak at the 
origin indicates the formation of chains that slide along the di¬ 
rector direction, thus preventing the development of strong po¬ 
sitional correlations between them. The morphologies found 
within the nematic phase in the system studied here are similar 
to the one observed for the smaller DSS (see section 3.2.1.). 



Fig. 7 Representative correlation functions for a GB-DSS mixture 
with G* = 1.5 andxr = 0.8: (a) two-dimensional pair correlation 

function in the uniaxial nematic (Nu) state at 

[(r*,p*) = (1.6,0.387)]. Distribution functions (b) (r^), (c) 

g|l^l (ry) and (d) g||^l ^ry^ for p* = 0.387 and various T*. 


At lower temperatures, the system forms a layered structure 
with highly interdigitated rods (see the snapshots in Fig. [^), 
also refiected by the light modulations of gy^l^ (ry) (see the 
black curve in Fig.[^). This is in contrast to the system with 
smaller spheres (cr* = 1.0) in which well defined layers are 
formed. Consequently, an increment of the diameter of the 
DSS strongly disturbs the translational order. An important 
finding is that the magnetic chains penetrate the smectic lay¬ 
ers rather than lying between the layers of rods (as it occurs in 


the corresponding GB-SS mixture forming an uniaxial smec¬ 
tic phase). The self-assembly of the DSS into magnetic chains 
alters dramatically the translational distribution of the DSS 
within the liquid crystalline matrix. It should be noted that 
the lamellar order is destroyed even when a lamellar config¬ 
uration for the GB-SS mixture is being used as initial con¬ 
figuration. Hence, a parallel arrangement of chains and rods 
is preferred instead of a perpendicular one. Furthermore, the 
inplane translational order of the rods indicates an isotropic 
positional arrangement characteristic of a SmA phase. 

3.2.3 GB-DSS mixture with a* = 1.7 We have further 
increased the diameter of the DSS particles to cr* = 1.7 and 
have studied the morphologies of systems consisting of A = 
720 particles at Vr = 0.8. Initially, we performed a cooling 
series starting from an isotropic phase. The system under¬ 
goes phase transformations to an uniaxial nematic state and, 
at lower temperatures, to a ’’glassy” state that does not pos¬ 
sess translational long-range order. This is a firm indica¬ 
tion that the DSS particles exert here a stronger perturbation 
on the LC matrix than for systems with smaller DSS sizes 
(7* < 1.5. It should also be noted that the volume fraction 
ratio 0r/0s — 2.44 is smaller than in the systems studied in 
previous subsections. 



Fig. 8 (a) State diagram of a GB-DSS mixture with cr* = 1.7 
(xr = 0.8), involving isotropic (I), nematic (Nu) and biaxial lamellar 
(Lb). The symbols indicate state points where the actual simulations 
were performed, (b) Global order parameters as functions of T* at 
p* =0.368. 


To get a more detailled insight into the behavior of the sys¬ 
tem we have also performed a heating series using lamellar 
structures obtained from the corresponding GB-SS mixture as 
initial configurations. A (r*,p*) diagram obtained for two 
different densities is shown in Fig[^. For clarity we describe 
the system at p* = 0.368. Interestingly, the lamellar structure 
persists also for the GB-DSS mixture. A visual inspection of 
representative snapshots (see Fig[^-b) reveals a microsepara¬ 
tion into alternating rod-rich and DSS-rich regions. Further¬ 
more, both species exhibit significant orientational order as re¬ 
fiected by increased order parameter values >0.7 (see Fig 
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8 3). The modulations of the (ry) function clearly indicate 

me formation of layers (see Fig[^-d). A pronounced interdig- 
itation of the rod species is also seen from the relatively slight 
modulation of g 11 ^1^ (r||) inFigj^. TheDSS self-assemble into 
ferromagnetic chains. As these are arranged in an antiparallel 
manner no net magnetization is exhibited. The most striking 
result is that the director of the rods is perpendicular to the 
director of the DSS with B We term this state biaxial 
lamellar Lb. 



Fig. 9 Representative simulation snapshots for a GB-DSS mixture 
with G* = 1.7, for Xr = 0.8 in the biaxial lamellar (Lb) state at 
[(L*,p*) = (1.3,0.387)] (a) both species are shown and (b) DSS 
particles are shown for clarity. The directions hr and hs are indicated 
by thick lines. Distribution functions: (c) g||^l and (d) 

g|l^l for p* = 0.387 and various L*. 


Upon further heating, the system undergoes a transforma¬ 
tion from a lamellar phase to an uniaxial nematic phase with 
respect to the rod species. Similar to the system with a* = 1.5 
there is a notable delay of the magnitude of the 5*^^^ order pa¬ 
rameter in comparison to 5*^^^ at the transition to the nematic 
state (see Fig[^). In summary, the above results indicate that 
(i) larger DSS particles destabilize conventional smectic or¬ 
der (either SmA or SmB) and (ii) the perturbation of the DSS 
particles to the LC host is such that the LC matrix does not in¬ 
duce, within the Nu state, any significant alignment to the fer¬ 
romagnetic chains either perpendicular or parallel to the mean 
orientation of the rods. 

3.2.4 GB-DSS mixture with cr* = 2.0 Here, we describe 
briefly the structure of a GB-DSS mixture with relatively large 


diameter, a* =2.0. For a detailed consideration of this system 
see Ref. |34l The self-organization of the magnetic particles in 
the nematic state differs in comparison to the mixture with 
(7* < 1.7. Remarkably, by decreasing the temperature starting 
from an isotropic state these systems exhibit spontaneously 
biaxial nematic ordering in which, on average, the director 
hr of the rods is perpendicular to the director hs of the DSS 
particles. Upon further cooling, a biaxial lamellar Lb is found 
similar to the GB-DSS mixture with a* = 1.7. In the Lb-state 
the ferromagnetic chains align, on average, perpendicularly to 
the director hr, and are arranged into an antiparallel manner 
resulting no net magnetization. 


3.3 Response to external magnetic fields 

3.3.1 GB-DSS mixture subject to external field with 

CT* = 1.0. 



Fig. 10 (a) Nematic Order parameter and (b) polar order parameter 
as a function of the external magnetic field strength for a GB-DSS 
mixture with G* = 1.0. Open symbols correspond to rods and solid 
symbols to dipolar spheres. 


In this section we examine the response of a GB-DSS mix¬ 
ture to an homogeneous magnetic field given by H = Hz 
at various field strengths W = mH/ksT. The field H is 
coupled to the permanent dipole m^ of particle i through 
the potential Ui = — ih^- • H. Three representative state 
points [(r*,p*) = (1.1,0.34)], [(L*,p*) = (2.4,0.40)] and 
[(r*,p*) = (1.2,0.40)] have been examined; the first two cor¬ 
respond to the isotropic and the third to the nematic phase of 
the field-free GB-DSS mixture (see the boundary lines in Fig 
The nematic and polar order parameters are shown in 
Fig as a function of the field strength. The director lis of 
the DSS particles is, on average, parallel to H. Therefore, the 
polar order parameter defined in Eq. ( p^ monitors the total 
magnetization of the system (i.e. the sum over the normal¬ 
ized magnetic moments). When the field is off in the isotropic 
state, for [(r*,p*) = (1.1,0.34)] and [(L*,p*) = (2.4,0.40)], 
the polar and the nematic order parameters are nearly zero. 
The nematic order parameter values for both the rod and the 
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DSS species increase by increasing the strength of the mag¬ 
netic field until saturation is reached at field strength ^ 15. 
Therefore, the rods species exhibit a field induced I-Nu trans¬ 
formation. The DSS particles form ferromagnetic chains that 
are oriented parallel to the direction of the magnetic field giv¬ 
ing rise to a net polarity that also saturates for even small 
magnetic fields (see Fig. [^. Remarkably, the linear ferro¬ 
magnetic chains induce orientational order to the rodlike par¬ 
ticles even for the isotropic state. The calculated nematic di¬ 
rectors of the species are on average parallel to each other and 
an uniaxial nematic phase is exhibited. Notably, the reverse 
phenomenon (in comparison to the systems described in sec¬ 
tion 3.2) is obtained here which has already been observed 
for real colloidal suspensions of rodlike and magnetic par- 
ticlesESl Under an external homogeneous field the magnetic 
spheres self assemble and form rodlike entities that induce ne¬ 
matic order to colloidal rods even in the isotropic phase. The 
most interesting finding is that the saturation value of the ne¬ 
matic order parameter of the rod species is larger for higher 
densities (p* = 0.40) than for lower densities (p* = 0.34). 
This means that the optical properties (such as birefringence 
in real systems) are expected to depend sensitively on the 
volume fraction of rods. Preliminary experimental resultsl^D 
for colloidal suspensions of rodlike particles with magnetic 
spheres indicate that the saturation value of birefringence in¬ 
creases by an increment of the volume fraction of rods at con¬ 
stant volume fraction of spheres. Finally, for the nematic state 
[(r*,p*) = (1.2,0.40)] the nematic order parameter is slightly 
increased under an external field (see Fig[T^). The polarity, 
on the other hand, increases significantly by changing from 
nearly zero to over 0.9 (see Fig[T^). 

3.3.2 GB-DSS mixture subject to external field with 
a* = 2.0. 



Fig. 11 (a) Nematic Order parameter and (b) polar order parameter 
as a function of the external magnetic field strength for a GB-DSS 
mixture with a* = 2.0. We use the parameter K defined as 
K = I — T*/T^ to indicate how far the system is from the I-Nu 
transformation temperature. For further details see the manuscript. 
Open symbols correspond to rods and solid symbols to spheres. 



Fig. 12 Representative simulation snapshots of a GB-DSS mixture 
with CTs* =2.0 (a) in the isotropic state for TT = — 0.9 and 
[(r*,p*) = (2.4,0.338)] and (b) same state point subject to an 
external homogeneous magnetic field of strength H* = 12.5. The 
direction of the field and the directors fir and fis are sketched as 
arrows. 


We have also considered systems with larger magnetic 
spheres (cr* =2.0) subject to an homogeneous magnetic field. 
More specifically, we have examined state points within the 
isotropic phase (at T* = 2.6 and T* = 2.4) and nematic phase 
(at r* = 2.2) at constant density p* = 0.338 [see (r*,p*) di¬ 
agram in Fig[T]i]. The nematic and polar order parameters as 
a function of the field strength are given in Fig [U we use 
the parameter K that is defined as ^ = 1 — to indicate 

how far away the system temperature is from the temperature 
where the I-Nu transformation occurs. 

Our results indicate (i) that nematic ordering is induced to 
the rodlike particles and (ii) the saturation values of the or¬ 
der parameters increase by decreasing the temperature of the 
system at constant density and (iii) the polar order parameters 
increase profoundly from zero (field off) and saturate to large 
values (over 0.8) even for small magnetic fields. The field- 
dependence of the order parameter for the DSS, on the other 
hand, is less affected by the temperature reaching similar sat¬ 
uration values at strong fields. 

One important finding is that the rodlike particles tend to be 
oriented perpendicular to the ferromagnetic chains (see Fig. 
[Tg. Hence, a biaxial nematic order can be induced within an 
initially isotropic state subject to an homogeneous magnetic 
field. 
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4 Conclusions 

In summary, using Monte-Carlo simulations we have stud¬ 
ied the influence of dipolar interactions and external magnetic 
flelds on the self-organization of binary mixtures of rods and 
soft spheres. We have implemented a tractable model of rods 
of Gay-Berne type (GB) and dipolar soft spheres (DSS) to ex¬ 
plore these types of systems. Due to the simple nature of the 
modelled interactions and particle shapes such prototypes of¬ 
fer conclusive stereotypes for addressing fundamental issues 
regarding the structure of complex liquid crystalline ferroflu- 
ids. 

In the first part of the paper we investigated the underlying 
reference system, that is, a binary mixture of rods and purely 
repulsive soft spheres (SS). For this GB-SS mixture we obtain 
fully miscible isotropic, nematic, smectic (SmB) and lamellar 
phases. A central result of the present study is that the ratio 
(7s / (To (diameter of the spheres compared to the width of the 
rods) determines whether the system undergoes a disordered- 
ordered phase transition to the smectic (cTs/cTo ~ 1 ) or lamel¬ 
lar phase (cT^/cTo > 1.5). In the smectic phase the spherical 
particles organize into linear arrangements along the nematic 
director. Moreover, we And that the spheres are fitting in the 
smectic layer of rods where the side-by-side rod-sphere ar¬ 
rangement is preferred. By increasing the size of the spheres 
the smectic order is destroyed in favour of a lamellar phase. 

The second part of the paper has been devoted to the impact 
of additional magnetic interactions on ordered states. To this 
end, we attached permanent magnetic moments to the centres 
of the soft spheres inducing strongly anisotropic long-ranged 
dipole-dipole interactions. For every configuration considered 
we observe a self-assembly of the dipolar spheres into chains 
due to a relatively large dipolar coupling (A > 3.5). The LC 
host stabilizes the orientational order of ferromagnetic chains. 
The ferromagnetic chains are randomly oriented up or down 
and the fluid has no macroscopic magnetization. Remarkably, 
the orientation of the dipolar spheres’ director with respect to 
the nematic director of the rods depends sensitively on the in¬ 
terspecies size-to-width ratio. For (Ts/cJo ^ 1.5 the ferromag¬ 
netic chains penetrate through the layers composed of (soft) 
rods forming an uniaxial smectic phase or at higher temper¬ 
atures an uniaxial nematic phase (parallel directors) and for 
ratios 1.7 the mixture forms a biaxial ferromagnetic 

lamellar phase (perpendicular directors). We note that for 
(7s/(To ^ 2.0 we have also observed a biaxial nematic phase 
at intermediate temperatures (for details see Ref. l34l) . 

As a further step, we have exposed the GB-DSS mixture 
to a (static) homogeneous magnetic held concentrating on 
isotropic states. Here we observe the reversed phenomenon, 
that is, the magnetic held induces orientational order to the 
chains (parallel to the held) that, in turn, imposes an orienta¬ 
tional order to the rods. Depending on the size of the dipolar 


spheres the rods are oriented (on average) parallel ((7s/(7o ^ 1) 
or perpendicular {g^/gq 2) to the held. The external held 
can thus be used to induce uniaxial or biaxial ferromagnetic 
order to a GB-DSS mixture which is isotropic in the absence 
of the held. 

There are several directions which we believe require fur¬ 
ther investigation. In the present paper we considered only 
relatively small spheres ranging up to diameters double the 
width of the rods. What would happen in the presence of even 
smaller spheres? Besides an increased tendency of demix¬ 
ing into rod-rich and sphere-rich regions these systems could 
possibly exhibit new morphologies. Another interesting point 
concerns the self-assembly considered here with further types 
of external influences, e.g., time-dependent magnetic flelds. In 
fact, a recent e^eriment on magnetic colloids under a rotating 
magnetic fleld^Sl revealed a phase transition from a disordered 
state to a layered state with hexagonal order; a phenomenon 
which was analyzed earlier in our groupie. For the present 
system we expect that this fleld-induced structural transition 
imposes a nematic order where the rods are aligned between 
the magnetic layers. This could be a promising route for de¬ 
signing functional liquid crystalline ferrofluids with potential 
magneto-optical applications. We therefore hope that our re¬ 
sults will stimulate further experimental investigations. 
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